OITS 668 
March 1999 



> 



CRITICAL BEHAVIOR OF HADRONIC FLUCTUATIONS 
AND THE EFFECT OF FINAL-STATE RANDOMIZATION 

! Rudolph C. HwaQ and Yuanfang Wu0 

^ I Institute of Theoretical Science and Department of Physics 

,—1 ; University of Oregon, Eugene, OR 97403-5203 

u • 

Abstract 

<■ 

^ • The critical behaviors of quark-hadron phase transition are explored by use of the 

Ising model adapted for hadron production. Various measures involving the fluctua- 
tions of the produced hadrons in bins of various sizes are examined with the aim of 
quantifying the clustering properties that are universal features of all critical phenom- 
ena. Some of the measures involve wavelet analysis. Two of the measures are found to 
I exhibit the canonical power-law behavior near the critical temperature. The effect of 

Q ' final-state randomization is studied by requiring the produced particles to take random 

0> ■ walks in the transverse plane. It is demonstrated that for the measures considered the 

i dependence on the randomization process is weak. Since temperature is not a directly 

I measurable variable, the average hadronic density of a portion of each event is used as 

^p^' the control variable that is measurable. The event-to-event fluctuations are taken into 

account in the study of the dependence of the chosen measures on that control variable. 
(— I I Phenomenologically verifiable critical behaviors are found and are proposed for use as 

• a signature of quark-hadron phase transition in relativistic heavy- ion collisions. 

X: 

1 Introduction 

If indeed quark-gluon plasma is to be produced in relativistic heavy-ion collisions at high 
energy, such as at BNL-RHIC, there must be a quark-hadron phase transition (PT) when the 
plasma becomes cool enough, at which point hadrons are created and subsequently detected. 
The properties of that PT have not been extensively investigated for a number of reasons. 
Theoretically, the physics is not perturbative, and the fluctuation of hadronic observables 
is beyond the scope of feasibility of lattice gauge calculations at present. Experimentally, 
there is the preoccupation with the conventional main-stream issues, such as the primordial 
signatures, hydrodynamical flow, interferometry, etc. The most serious obstruction, however, 
is the common belief, based mostly on an extrapolation from low-energy nuclear collisions, 
that a hadron gas surrounds the plasma and that the final-state interaction in that gas 
obliterates all features of the PT, leaving the hadrons with no significant memory of their 
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birth to be registered at the detector. One of the aims of this paper is to investigate the 
type of hadronic signature that can survive the final-state interaction and thereby to dispel 
the myth that critical behavior of quark-hadron PT is not observable. 

If we take statistical physics as a guide, then we should expect the critical phenomena 
related to quark-hadron PT to be a vibrant area of study. The familiar critical behaviors 
of conventional systems are, for example, C ~ |T — T^p" for the heat capacity of a thermal 
system, or x ~ \T — Tc\~'^ for the susceptibility of a magnetic system, where is the critical 
temperature. Is there any similar behavior for a quark system in its transition to hadrons? 
Can we find an observable in heavy-ion collisions, call it /i, such that its behavior near the 
critical point is 



Will /i survive the randomization process in the hadron gas, if it exists? Since the tempera- 
ture T is not directly observable in heavy-ion collisions, is the critical exponent k measurable? 
If not, what is the nature of the critical behavior that is subject to experimental test? These 
are some of the questions that we shall address and find answers to. 

Our investigation cannot be based on nonperturbative QCD, which is too difficult to 
implement in treating the problem posed. In the absence of a reliable dynamical framework 
that can substitute for QCD, we appeal to the lesson already learned in statistical physics 
in that the universal properties of the critical phenomena are independent of the details 
of the microscopic dynamics. Only dimension and symmetry are of most importance. In 
the Ginzburg-Landau theory a universal scaling behavior on multiplicity fluctuations was 
derived and later verified experimentally in quantum optics at the threshold of lasing 
@. That being a mean-field analysis, a more elaborate investigation that involves the study 
of spatial fiuctuations near the critical point was subsequently carried out in the framework 
of the Ising model P]. The result confirms the average behavior found in [U, but the work 
goes beyond the mean-field approximation by examining the PT process on a 2D lattice in 
a way that can best simulate the hadronization process. In this paper we continue to adopt 
the same approach, the only one available short of QCD, to study the fiuctuation of spatial 
patterns of hadronization during PT, and to investigate whether the measures chosen can 
survive the randomization process of final-state interaction. 

Instead of working with factorial moments, which are the measures of choice in the 
study of intermittency we shall consider elementary wavelet analysis to enhance the 
effectiveness of describing clusters of all sizes. Furthermore, we shall go beyond intermittency 
and study the fluctuations of spatial patterns [|, 01 • A measure that is based on the wavelet 
coefficients will be investigated in our search for a description of the critical behavior. 

Our physical picture of the heavy-ion collision process at high enough energy and nuclear 
masses for quark-gluon plasma to be formed is the conventional one. After formation, a 
cylinder of high temperature plasma that is created undergoes rapid longitudinal expansion 
and slower transverse expansion. The interior is at T > Tc so the medium stays in the 
deconfined state during most part of the time evolution of the system. We assume that 
with two fiavors the quark-hadron PT is second order and that it takes place on the surface 
of the cylinder where T is lowest, namely or slightly below. Thus the geometry of the 
system undergoing PT is 2- dimensional. For a time interval of the order of 1 fm/c, which 
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is roughly the hadronization time for the formation of a hadron, we map a section of that 
cyhnderical surface to a square lattice to be studied in the Ising model. Hadrons are not 
formed uniformly on the surface. At different time segments of the evolutionary history the 
spatial patterns of hadron formation are different. Their fluctuations are simulated by the 
changing configurations of the Ising lattice. It is on the basis of this mapping between the 
physical surface and the mathematical lattice that we pursue our extensive investigation of 
the 2D Ising model in this paper. 

2 Hadron Density in the Ising Model 

Our first task is to establish a quantitative connection between the hadrons produced in a 
heavy-ion collision and the spins on a 2D lattice in the Ising model, which we use to simulate 
a second-order PT. A number of intermediate steps are needed to make that connection; 
they have been discussed in Ref. [0], which we summarize here. 

The basis of our approach is the Ginzburg-Landau (GL) theory of PT. On the one hand, 
there is the Ising system of spins that can implement the GL description through a simple 
dynamical model capable of generating microscopic fluctuations 0. On the other hand, 
particle production in the form of photons is well described by the GL theory, when a laser 
is operated at the threshold of lasing ^ . Our prediction for the scaling behavior associated 
with hadron production in the GL formalism [|l], |lOl turns out to be verified experimentally 
by photon production in quantum optics Thus we shall assume that the relationship 
between the GL order parameter and the hadron density to be the same as that for the 
photon density in the laser problem which is rigorously known. 

The GL free energy density is 

T[(t>\=a\ct){z)\^ + h\ct){z)\'' + c\dct)/dz\' , (2) 

where (f){z) is the order parameter for the spatial coordinate at z in a 2D space. The 
parameters a, h and c need not concern us here. (I){z) is the link between hadrons and Ising 
spins. On one hand, the hadron density, in analogy to the photon density, is 

p{z)^\m? , (3) 

while on the other hand (j){z) is the local mean magnetization field, i.e., 

m = K"'T.s^ > (4) 

where the sum is over all lattice sites in a cell of area = e^. In the absence of external 
field the total magnetization of the whole lattice may be either positive or negative. We 
define Sj to be positive if the spin is aligned along mi- Thus if aj is the lattice spin having 
values ±1, so that = J^jeL^ then we define 

Sj = sgn (mi) aj , (5) 

where sgn (m^) stands for the sign of m^. This change of sign when changes sign is 
irrelevant to the connection between the GL theory and the Ising model; however, is 
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important in establishing a connection between hadron multiplicity and Ising spin. When 
T > Tc, hadronization is inhibited, and (j){z) is near zero since the summation of spins in 
a cell in (^) tends to cancel the site spins in the disordered state. If it is nonzero due to 
fluctuations, hadrons can be formed, but we associate them only with </)(z) > 0, as we would 
for the ordered state, T < Tc, where hadrons are formed more copiously. In order to relate 
different event of collision with different configurations of the Ising spins in a uniform way, we 
must associate hadron formation with spins aligned with m^, i.e., J2j Sj > 0, not J2j > 0. 
To summarize, the hadron density in the ith cell is 



j€Ac(i) 



E^. ' (6) 



where B is the Heavyside function. The factor A in (|^) relates the mathematical quantity on 
the right side defined on an Ising lattice to the physical quantity of the number of particles in 
a cell emitted from the surface of quark-gluon plasma formed in a heavy-ion collision. Since 
A is an unknown parameter in our attempt to simulate quark-hadron PT by the Ising model, 
we shall vary A and regard any results that are sensitively dependent on A as unphysical. 

In our Monte Carlo simulation of the Ising model we use the Wolff algorithm p|, |Tl| to 
calculate the spin configurations for a lattice of size L = 256, using the Hamiltonian 

H = -jJ2a,a, , (7) 

where the sum is over the nearest neighbors. We choose the size of a cell to be e = 4 so 
the maximum value of p/A in the ordered phase is (4^)^ = 256. Thus A would have to be 
quite small to render the multiplicity in a cell to be a reasonably small number. Near Tc we 
expect p to vanish in many cells. A bin of size 6 consists of Nc = cells, and the whole 

lattice consists of M = bins. The average density, averaged over all cells bins /c, 

and configurations e, is 

(p) = j^Y.MY.jfY.p^ ' (8) 

J\l IVl l\c j=i 

where M is the number of simulated configurations. In the following the results obtained 
are for M = 10^. 

The dependence of (p) on T is shown in Fig. 1. It is evident that a continuous PT occurs 
at Tc somewhere between 2.0 and 2.5 in units of J/ks- The precise value of Tc has been 
determined in ^ by examining the scaling behavior of the normalized factorial moments Fq. 
It is found that Fg behaves as M'^'' only at T = 2.315. Since the critical point is characterized 
by the formation of clusters of all sizes in a scale independent way, we identify the critical 
temperature at Tc = 2.315. Note that (p) does not vanish for T > Tc, contrary to the 
behavior of the average magnetization in the limit of vanishing external field. The difference 
is due to our definition of hadron density p in (j^) as opposed to the average magnetization 
being (m^). For heavy- ion collisions it is reasonable to expect that hadrons can be produced 
even at T > Tc, though at reduced rate, due to local fluctuations in quark density and 
interactions. 
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For the sake of visualization of the spatial pattern of hadrons formed at Tc, we show 
in Fig. 2 a particular configuration on the lattice simulated at Tc. Squares of various sizes 
denote proportionally the cell densities of hadrons. It is evident that clusters of all sizes are 
formed. That is the well-known phenomenon at PT. Here it is the hadron density in a cell 
that is shown, not just spin up or down in the Ising problem. 

Now, in search for a behavior analogous to (|l]), we focus on the T dependence of (p) for 
T < Tc. In Fig. 3 we show (p) — (pc) vs — T in a log- log plot with A set at 1, (pc) being 
(p) at Tc. For the range of T examined, it is clear that there exists a power-law behavior 

(p) - (pc) oc (Tc - T)" (9) 

with 

ri = 1.67 . (10) 

For comparison, the susceptibility of a magnetic system has a critical exponent of ~ 1.3. 

If the final-state interaction in the hadron gas exists for several fm/c, long compared to 
the typical time for resonance decay (~ 1 fm/c), the multiplicity of stable hadrons (in strong 
interaction) produced would remain unchanged during that phase, since a temperature of 
~ 140 MeV is not high enough to cause further production of particles. Thus the average 
hadron density is insensitive to the final-state interaction. However, (y) is not what is usually 
regarded as critical behavior, for which the exponent in a formula such as (|ID is negative. 
Our search for a suitable measure of critical behavior therefore continues. 



3 Fluctuations of Hadron Density 

Since the fluctuation in particle production is large near the critical point, as evidenced by 
Fig. 2, we next consider the fluctuation of p from the average (p). Let the average density 
in bin k of size 6 be 

p^=y (11) 

where the sum is over all cells in bin k, the total number of which is = {5/ef'. Near Tc, 
Pfc fluctuates from bin to bin, especially for small 5. Let us consider the moments 

where (p) is given in (^. By studying the ratio pkj (p), we avoid the dependence on the 
unknown A; moreover, we have Aq = 1 aX q = 1. A representation of the properties of Aq is 
given by 



J(M) = —Aq 



d_ 

dq 



q=l 
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One can also study q = 2 and higher powers, but J{M) is sufficient for our purpose here. 

In Fig. 4 we show the result of our simulation on J(M). If Aq had a power- law dependence 
on M = (L/5) , then J{M) would depend linearly on In M. Evidently, Fig. 4 does not exhibit 
any linear dependence at any T. However, the behaviors at various T < Tc are similar. The 
situation is reminiscent of the behavior of Fq in the GL theory where strict scaling 
Fq oc M'^'' is not valid, but Fq oc is valid. Similarly, we can consider here the dependence 
of J{M) on Jc{M), where Jc{M) = J{M,T = T^). As can be seen from Fig. 5, the points 
for T < Tc can be well ffited by straight lines. Let us therefore write 

J{M,T) = a{T)JciM) . (14) 

This is the property of a generalized form of scaling [|l^, |l^ 

Aq{M,T) ^ giMY"^^') , (15) 



where g{M) is some function of M independent of T. It follows directly from ( |13|) that 



aiT) = r'iT)/r'iTc) , (16) 



where t'(T) = dTq{T) / dq\^^^. We do not pursue the study of the generalized scaling (|15D 
here. 

From the slopes of the straightline fits in Fig. 5, we can exhibit the T dependence of 
a(T), shown in Fig. 6. We display a{T) in log-log plot so as to explore the possibility of the 
critical behavior 

a(T)oc(T,-T)-^ . (17) 

We have no illusion that ( [T7| ) could be valid in the limit T T^, since a{Tc) = 1, by 
definition. However, for a range of T below Tc, Fig. 6 does show an approximate linearity. 
The corresponding value of the critical exponent rj is 

C = 1.88 . (18) 

Since (|3) is invalid in the immediate neighborhood of Tc, it cannot be taken seriously as the 
defining characteristics of quark-hadron PT. Instead of (p!?]), we regard Fig. 6 as containing 
more information about the fluctuations of the hadron density near the critical point. How 
such a behavior can be checked experimentally will be discussed in a later section below. 



4 Fluctuations in Clustering 

In the previous section we have studied the fluctuation properties of hadron density from 
bin to bin for all events. Since the bin density pb is an average of the cell densities in a bin, 
it is not sensitive to the fluctuations of the cell densities within a bin. Clustering being such 
a characteristic feature of critical behavior where clusters of all sizes can occur, we need to 
go beyond ph to capture the fluctuations of multiplicities among neighboring cells within the 
same bin. To that end we shall work with multiplicities instead of densities of hadrons, not 
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only for calculating spatial fluctuations, but also for considering temporal fluctuations in a 
later section when randomization due to final-state interaction is taken up. 

In terms of the cell density pi defined in the cell multiplicity is simply e^pi. Since 
A in (ID is unspecified, let us absorb the factor into A, i.e. redefining e^A A, and write 











rii = I 


A 















(19) 



where I[x\ signifies the largest integer < x. Thus is discrete, although A is treated as 
a continuous parameter < 1. For every configuration simulated on the lattice we have a 
matrix C (ci, C2) that gives the value of the cell multiplicity rii at the cell position (ci, C2) in 
the 2D space, where q ranges from 1 to £ = L/e; i.e., C (ci, C2) = ni[i = (ci, C2)]. Figure 2 is 
a pictorial representation of a possible configuration whose numerical description would be 
a particular matrix C (ci, C2). 

Now, for every bin in that £ x £ space we consider a transform through a function ip^ 
that has the following property. Divide a square bin into four quadrants: in the two diagonal 
quadrants = +1, in the two off-diagonal quadrants ip^ = — 1. If we use x and y to denote 
the two orthogonal coordinates of the square bin, normalized to unit length on both sides, 
then the precise definition of ip^{x,y) is 

,y) = l, 0<x<l/2, 0<y<l/2, 

= -1, l/2<a;<l, 0<y<l/2, 

= -1, < X < 1/2, l/2<y< 1, (20) 

= 1, l/2<x<l, l/2<y<l, 

= 0, elsewhere . 



This is essentially a 2D Haar wavelet ||^, but any prior knowledge about wavelet analysis is 
not needed for our discussion here. 

We label a specific bin in the ix i space of cells by a triplet (j, ki, k2), where the bin size 



6j is 



5, 



(21) 



and the location of the bin is specified by {ki, k2), with ki ranging from to 2-' — 1. Clearly, 
the total number of bins is 2^-' . Thus j gives the level of zooming, and (fci, k2) describes the 
shift at that level. The label j used here and in the following should not be confused with 
the symbol j for site position on the lattice in the equations from (4) to (7). We can now 
define our basic wavelet function in the i"^ space 



i'jMM (ci, C2) = yj- - ki,j- - k2 



(22) 



It vanishes everywhere except in the jth level bin located at (fci, ^2), and in that bin it can 
only be ±1 depending on the quadrant that Ci and C2 fall into. 
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Using ipj,ki,k2 (ci, C2) we define the wavelet transform of a configuration of cell multiplici- 
ties 



'^jMM = i'jMM (ci, C2) C (ci, C2) , (23) 



Cl,C2 



where the sum is over all cells. It is clear that if C (ci, C2) is uniform in a bin (j, ki, /C2), then 
the wavelet coefficient for that bin, Wj^kiM^ zero, although the average density ph for that 
bin is not. Thus Wj^kiM measures the fluctuation from the average in a specific way. By 
varying j we can learn a great deal about the landscape of C (ci, C2). 

To study the fiuctuations of the clustering within each configuration, but averaging over 
all configurations, we investigate the scaling properties of the moments of the wavelet coef- 
ficients, Wj^kiM- "^^^ wavelet coefficients, being more sensitive to fiuctuations than Fg, have 
a better chance of generating a discernible signature of clustering patterns formed near the 
critical point. To that end, we define for a fixed j, for which M = 2^-^, 

B,{M) = ((^^y), (24) 



{w) 



where 



I -1^ I 2^-1 



e=l ki2=0 



Furthermore, we have 



K{M) = —Bg 



d_ 

dq 



9=1 



^^In^^V (26) 



which is the quantity that we shall investigate. In studying the ratio Wj^^iMl ('^)' 
effects of A in (|19]) tend to cancel, but not completely because rii is a discretized version 
of a noninteger quantity at each cell. Thus some dependence of K{M) on A is expected, a 
property that will be examined. 

The erraticity analysis originally designed to describe the event-to-event fiuctuations 
was done using factorial moments Fg. For the problem of PT it was found that the erraticity 
indices turn out to be very small (of order 10^'^) |]T^, smaller than what experimental data can 
probably resolve. What is proposed in (^)-(^) is not erraticity analysis, and is not aimed 
at studying the fiuctuations from configuration-to-configuration. Instead, Bg{M) and K{M) 
are measures of the fiuctuations from bin to bin that are averaged over all configurations. 

In Fig. 7 we show K{M) vs In M for a number of values of T < T^. The dependence on 
M is evidently far more severe than that of J{M) at large M. Due to the fact that j must 
be an integer, we cannot study at smaller steps in M. The fiuctuations in small bins are 
clearly very strong, and they do not vanish in large bins, lest K{M) would approach zero, 
as M — > 1. It means that hadrons are formed in clusters of all sizes, roughly with as much 
empty space on the lattice (unhadronized quarks in real space) between the clusters as there 
are hadronized cells. 
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To facilitate a description of the T dependence of K{M, T), we plot K{M, T) vs Kc{M) = 
K{M,Tc), varying M parametrically, as shown in Fig. 8. Despite our inability to add points 
corresponding to intermediate bin sizes, the linearity of the behavior is self-evident. The 
points can be well fitted by straight lines, yielding the form 

K{M,T) = f]o{T) + f3{T)K,{M) . (27) 

Thus we can extract the slope /3(T) that depends only on T, and not on the bin sizes. It 
therefore characterizes the intrinsic properties of the system under PT, and may be more 
general than the Ising model used to derive it. 

Figures 7 and 8 show the results of the calculation where A has been set to 1 in ([T9|) . We 
now decrease A to see how P(T) depends on A. That is shown in Fig. 9. Considering the fact 
that A has been varied over two orders of magnitude, the form of (3{T) is basically stable. If 
we fit the last three points by straight lines in accordance to the formula 

/3(T)oc(Te-T)-" , (28) 

the value of k varies weakly with A, as shown in Fig. 10, except at A = 0.02 and 0.03 for 
unimportant reasons having to do with the discretization of rii. Hereafter we shall adopt 
A = 0.05 as the representative value and set 

K = 2.2 . (29) 



Although (28) is not valid as T — Tc, since by definition (3 (Tc) = 1, its approximate 
validity in a narrow range of T nevertheless captures a behavior remarkably similar to the 
conventional critical behaviors of familiar thermal and magnetic observables. 

Since a{T) and /3(T) look very similar in Figs. 6 and 9, it is sensible to relate them 
directly, using T only as an implicit variable. That is, we plot (3 against a, as in Fig. 11, 
where all of the five temperature points are shown. They correspond to the A = 0.05 points 
of a(T) and P(T) in Figs. 6 and 9. We see that there exists a linear relationship, as exhibited 
by the straightline fit 

(3 = a + ba , (30) 

where a = —0.0623 and b = 1.11. It is particularly remarkable in view of the nonlinear 
behavior in Figs. 6 and 9. There has been no prior suspicion that such a relationship should 
exist, so the experimental verification of (^) would be of great interest. 

One may wonder how (|3y) can be consistent with (^^ and (^8]) that have different 
exponents. Indeed, if one examines the linear fits in the log-log plots in Figs. 6 and 9, the 
deviations from the straightlines are larger than the deviations of the points in Fig. 11 from 
the solid line in that figure. In other words, when discrepencies from best fits are taken into 
account, the linear behavior exhibited in Fig. 11 is a better fit of the dependence of /3 on a, 
which is especially significant in view of the full range of all five points being related. 

In the discussion in connection with ( PT| ) and ( pSj ) we have emphasized the behaviors 
of a and (3 in the region where T is not in the immediate neighborhood of T^. We now 
examine further the region where T is closer to Tc. In Fig. 12 we have plotted the curves 
for K{M, T) vs J(M, T) for T < Tc. The solid curve for T = Tc is of particular significance 
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if the PT occurs only at Tc- It is evident that the curves for T close to Tc have very 
similar shapes. Thus it appears that by shifting those curves they can all be overlaid on a 
universal curve that coincides with the curve at Tc. Since shifting horizontally and vertically 
correspond to rescaling the horizontal and vertical variables, we have replotted the curves 
in a set of new variables that involve changes in scales. The result is shown in Fig. 13, 
where K{T){T/Tc)'^ is plotted against J{T)(T/Tc)^^^ , with M being varied implicitly. All 
the curves for 2.28 < T < T^. = 2.315 overlap and form a "universal" curve. Such a scaling 
behavior involving a range of T may be of theoretical interest only at this point, since T is 
not measurable. Nevertheless, the existence of such a universal curve is suggestive of some 
property that is like the Kadanoff scaling, where a change in the size of spin blocks can lead 
to a shift of the effective temperature [|l^. Of course, the quantities involved here are quite 
different, and the reason for the universality found here deserves further investigation. 



5 Observable Consequences 

The results that we have obtained so far from analyzing the data generated by the Ising model 
are the T dependences of (p), oi{T) and /3(T). Such results have no direct phenomenological 
implications, since the temperature T is not a directly measurable quantity in heavy-ion 
collisions, let alone the precise value of Tc. Thus equations such as (H), (|l3) and ( p8|) 
are primarily of theoretical interest. They are the results of a search for behaviors that 
resembles the conventional critical behaviors. Now we ask, given these theoretical hints, 
what observable consequences of PT can be suggested that are experimentally testable. 

To proceed, it is necessary to recognize first the important aspects in which the realistic 
situation of quark-hadron PT in heavy-ion collisions is more complicated than what has 
been simulated on the Ising lattice. A lattice configuration generated by the Ising model 
can represent the hadronization pattern on the surface of the plasma cylinder for only the 
time duration of the hadronization time th, roughly 1 fm/c. As time proceeds, the pattern 
changes, just as the lattice configuration changes upon repeated simulation. Under the 
assumption that the plasma interior is hot, with T > T^, phase transition occurs on the 
surface at T < Tc over a period of time long compared to th- What needs to be analyzed 
in the heavy-ion experiment is the hadronization patterns separated from one another in 
time segments of about th long. While that may be hard to achieve in present experiments, 
one can approximate that by making cuts in small Apx intervals, since one expects from 
hydrodynamical considerations that a correspondence between px and evolution time r exists 



such that hadrons emitted early have higher pt than those emitted later . We proceed 
under the supposition that it is possible to obtain from the heavy-ion data collected at a 
fixed transverse energy Et many samples consisting of portions of the events that correspond 
to the hadronic patterns in the ?7-0 space and are the realistic counterparts of the simulated 
hadronic configurations C (ci, C2) on the Ising lattice. 

In the Ising model we can adjust the temperature T and study the T dependences of 
various measures. In heavy-ion collisions that cannot be done. We then consider two possible 
scenarios. First, we assume that the quark-hadron PT occurs at Tc only, and no other T. 
Second, we consider the possibility that various conditions of hydrodynamical flow may lead 
to hadronization at various T below Tc at different points of time in the evolutionary history 
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of the plasma. Let us refer to these two scenarios as A and B, respectively. 

In scenario A we give up the possibility of studying the T dependence of PT. Our analysis 
of the data generated on the lattice nevertheless can serve as a guide to the analysis that 
can be done on the heavy-ion data. Since Jc(M) and Kc{M) are both observable, their 
dependences on M in Figs. 4 and 7 can both be checked. Due to the absence of an absolute 
relationship between bin sizes on the lattice and those in the experiments, it is best to 
examine the dependence of on with M being parametrically varied. Such a dependence 
is already shown by the solid line in Fig. 12 and by Fig. 13. Clearly, such a behavior is 
independent of the definitions of bin sizes on the lattice or in the experiment, and is a 
definitive characterization of the fluctuations of hadronic clusters at PT. The verification of 
the relationship bewteen Kc and Jc would imply that the hadrons detected are created by a 
quark-hadron PT at Tg. 

In scenario B where quark-hadron PT can occur at a range of T, the problem is signifi- 
cantly more complicated and will be the subject of our discussion for the remainder of this 
section. Let us first define a sample to be a portion of an event obtained by certain cuts 
(such as by a narrow Apx cut discussed in the beginning of this section) so that a sample 
corresponds to a configuration simulated on the Ising lattice. We assume that a sample is 
the result of hadronization at a common T, which may differ from T^. That T may vary 
from sample to sample. For every sample the average density p of the hadronic configuration 
can be determined experimentally and should correspond to the average density calculated 
on the lattice according to 



the latter is averaged over all configurations at the same T, whereas p refers to a specific 
configuration on the lattice or a specific sample in the experiment. Since p can vary from 
sample to sample (as they do from configuration to configuration on the lattice) even when 
all samples are at the same temperature, let us first examine the distribution P(p, T) of p at 
fixed T. We have calculated that distribution at three representive values of T with the result 
shown in Fig. 14. As T is decreased from Tc, the average density increases, as more hadrons 
are produced in a phase transition. The shapes of the distributions vary quite significantly 
as functions of T. On the lattice the cell densities are bounded by a maximum, according 
to (|^), of Pi,max = 256A = 12.8 (for A = 0.05). The normalization of p should not be taken 
seriously in phenomenology, since the density on a lattice cannot be rigorously related to the 
hadron density in an experiment. We should therefore regard p as having an arbitrary scale. 
Nonetheless, these distributions at fixed T are important evidences that the event-to-event 
(and therefore configuration-to-configuration) fluctuations are large even if all kinematical 
variables of the collision processes are fixed, and if the uncontrollable dynamical variable T 
at hadronization is held constant. 

In view of the rather wide distributions in Fig. 14 (which are all normalized to 1), it is 
clear that any experimental measure categorized according to the observed p must involve a 
convolution over the unobserved T. Thus, for example, in studying K{M) we should convert 




(31) 



P 



This p is not the same as the theoretical overall average density (p) defined 
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its T dependence to p dependence by use of the formula 

K{M, p)=JdT P{p, T) Q{T) K{M, T) , (32) 

where Q{T) is the probabihty that the quark-hadron PT will take place at T < Tc. 

We have no deep insight on what Q{T) should be. It can depend on the thermodynamics 
and hydrodynamics of the heavy-ion collision problem. One would expect that PT should 
take place in the immediate neighborhood of T^. We shall take an exponential form for 
Q{T) as a working hypothesis. Since the calculations in Sees. 3 and 4 are done for the range 
of T between Tc = 2.315 and T = 2.25 and interesting critical behaviors have been found 
theoretically in that range, we shall proceed on the assumption that Q{T) is important only 
in that range. Accordingly, we adopt the following form 

Q(T) =go exp[-(Tc-r)/0.08], 2.24 < T < 2.32, (33) 

and require that it is zero elsewhere. The normalization factor is go = [0.08(1 — e^^)]^^. 

We now can redo the calculations on J(M) and K{M) as in Sees. 3 and 4 by use of the 
convolution as in (|32D and present the results in terms of p. We shall set A = 0.05. In Fig. 
15 we show J(M, p) for various bins of values of p. Note that the highest curve (in solid 
line) is for 3.75 < p < 4.0, for which we shall use po = 3.88 to denote that bin, a value that 
shall play the role of in Fig. 3. At lower p the curve for J{M) is lower; this fringe effect 
shall be ignored in the scaling analysis below. That is, we shall consider only p > pq in the 
following. The similarity of the curves suggest following the same procedure as in Sec. 3. 
Thus we define Jo{M) = J{M, p = po), and examine J(M, p) vs Jq{M). The result is shown 
in Fig. 16. The linear behavior is remarkable. Defining the slope to be a(p), we show its p 
dependence in Fig. 17. A power-law behavior can be identified, as indicated by the straight 
line, yielding 

a(p)oc(p-po)-^", C = 0.97. (34) 

The same procedure can be applied to K{M,p), as done in Figs. 18-20. With Kq{M) = 
K{M, Po) the linearity of the lines in Fig. 19 is even better than that in Fig. 7. Again, 
defining the slope to be /3(p), we obtain from Fig. 20 the behavior 

/3(p)oc(p-po)"^ R = 0.97. (35) 

Evidently, the behavior of f3{p) is essentially identical to that of a{p). Equations (0) and 
(p5|) represent the closest that we can identify as the critical behavior of quark-hadron PT. 
This time all the quantities involved are experimentally measurable. 

Finally, we can plot (3{p) vs a(p) with p being the parametric variable. All five points 
fall on a straight line in Fig. 21. The linear behavior can be well expressed by the amazing 
formula 



= Pip) 



with unit slope. We have no explanation for the simplicity of (|36|) . 
significant if this behavior can be verified by experiment. 



(36) 

It would be highly 
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6 Effects of Randomization 



What we have studied so far are the properties of the hadrons formed on the cyhndrical 
surface, which is mapped onto the Ising lattice. The observed hadrons that reach the detector 
must traverse the space between the cyhndrical surface and the detector. If a hadron gas 
surrounds the plasma cylinder, then a hadron formed on the surface must undergo final- 
state interactions with other hadrons before it can move as a free particle to the detector. 
Assuming no further particle production in this last phase, the hadrons can only shift their 
positions in the //-^ plane in a random way, since the hadron gas plays the role of random 
scatterers and can offer no organized forces on a traversing hadron. 

We shall take into account this randomization process by requiring each hadron to take 
random walks on the lattice. The step size is the distance between neighboring cells, i.e., e. 
If a cell has > 1, then each particle in the cell takes steps that are independent of what 
other particles in the cell do. A particle has five possible positions in each random step, four 
neighboring cells plus the original cell for no transverse movement. Denoting the number of 
steps by z/, we consider v ranging from 1 to 6, although it is generally believed that 3-4 steps 
are adequate to represent the effect of randomization by the interactions in a hadron gas. 

In considering random walk we must focus on the whereabout of the particles, so hadron 
density must be converted to particle number by a specific choice of A. Although we have 
avoided that choice in Sec. 3, we have found in Sec. 4 that A = 0.05 is an appropriate choice. 
Thus we use the same A now in calculating the effect of random walk on J(M, T) and then 
on q;(T). Fig. 22(a) summarizes the effect for v up to 4, where we stop to avoid overlap of 
the points on the figure at higher v. The general shape of a(T) is unchanged. Quantitative 
dependence of the index defined in (|T7|), on v is shown in Fig. 22(b), exhibiting insensitivity 
to randomization. 

The same analysis has been done for K{M,T) and then (3{T) with the result shown in 
Fig. 23. In this case the effect of random walk is more pronounced, especially at low z/, but 
for 4 < u < 6 the effect stabilizes. Although (3{T) suffers more change than a(T), the overall 
effect does not invalidate the usefulness of the measures that we have considered. We show 
in Fig. 24 the effect of u steps of randomization on the dependence of (3 on a. Evidently, the 
general linear dependence is unaltered by the randomization process. 

We shall not redo the above calculation for a and j3 as functions of p, since the effect of 
randomization should be similar on the experimental a{p) and /3(p). If the experimental data 
can produce a figure such as Fig. 21, showing that (3{p) is essentially the same as a(p), we 
can conclude on the basis of the study considered here that the effect of randomization is not 
so severe as to render the whole approach ineffective. We infer that the measures J(M) and 
K{M) must not be very sensitive to the redistribution of particles in the transverse plane, 
since it is not the individual local property of a hadron but the scaling properties of clusters 
of hadrons over a wide range of bin sizes that the fluctuation measures have extracted. 

7 Summary 

Using the Ising model and defining hadrons on the lattice in accordance to the Ginsburg- 
Landau description of quark-hadron PT, we have found a number of interesting features 
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about the critical behavior. First of all, the average hadron density on the surface of the 
plasma cylinder for a time interval around the hadronization time of a hadron depends on 
the surface temperature according as 

{p)-{p,)^{T,-TY , r/ = 1.67 . (37) 

Then there are two quantities Q;(T) and l3{T), defined in (|T^) and (p?]), that behave as 

a(T) oc (Tc-T)"^ , C = l-88 . (38) 
/3(T)oc(T,-T)-^ , K = 2.2 . (39) 

These are behaviors of theoretical interest, since T is not directly measurable. Nonetheless, 
(^) and (|39|) confirm that there exist measures in the hadronic variables that exhibit the 
canonical form of critical behaviors. 

Since the average density p in a sample is measurable, we have shown how the critical 
behaviors can be reexpressed in terms of p as 

a{p) (X (p - po)-^ /?(p) (X (p - po)-", (40) 



( = R = 0.97 . (41) 

The equations in (^) are expressed in terms of observables and should be checked by ex- 
periments. Furthermore, the direct dependence of (3{p) on a{p) over a wider range of p can 
also be examined. If it turns out that they are related approximately by a(p) ~ /3(p), then 
that would be a less restrictive, yet nevertheless valid, experimental signature of the critical 
behavior of a quark-hadron PT. 

Experimental study of J{M) and K{M), defined in (|T3|) and (^), involves less analysis 
and thus are more reliable phenomenologically. If PT occurs at T = Tc only, then Jc{M) 
and Kc{M) are directly measurable. The verification of Kc vs Jc, as shown by the solid line 
in Fig. 12, would strongly indicate the occurrence of a PT. 

None of the phenomenology suggested above would mean anything realistic if the hadron 
gas in the final state erases all memory of the properties of phase transition before the 
particles reach the detector. We have studied the effect of the randomization due to hadron 
gas on our measures, and found that within the framework of our modeling none of them are 
seriously influenced by the randomization procedure so long as the number of random steps 
in the ?7-0 plane is not too large. Those measures involve the fluctuations from bin to bin of 
either the average density of hadrons or the wavelet coefficients of the spatial patterns, both 
of which are evidently insensitive to the randomization. 

An important assumption is made in this work, namely: cuts in small Ap-r bins can be 
made to select small At bins in real time that can exhibit the various hadronic configurations 
of a phase transition process. Based on the intuition gained in this work, we conjecture that 
the general features of our result are insensitive to the precise validity of this assumption, 
provided that ApT is small enough. Only by doing the appropriate analysis of the actual 
heavy-ion collision data can we assess the plausibility of this assumption and the feasibility 
of the whole program discussed here to extract the critical behavior of quark-hadron phase 
transition. 
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It should be noted finally that there are two aspects about the work presented here. 
One is about the theoretical understanding of the properties of the critical behavior at the 
quark-hadron PT. The other is the attempt to find the means to determine those properties 
experimentally. Even if the latter effort fails because of our incomplete understanding of 
the actual collision and hadronization process, the former will remain valid and stand ready 
for verification by better phenomenology to come. Our result on the insensitivity of our 
proposed measure to the final-state effects provides the encouragement to persist in this 
difficult but worthwhile quest. 
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Figure Captions 

Fig. 1 Temperature dependence of the average hadron density in arbitrary scale. 

Fig. 2 A representative configuration of hadron production in rj-cj) plane as simulated on 
the Ising lattice of size ixi, where £ = L/eis the number of cells along each dimension. 
The size of a square is proportional to the hadron density of a cell. 

Fig. 3 Temperature dependence of average density for T < T^. t] is the slope of the 
straightline fit. 

Fig. 4 The dependence of J{M) on M for various values of T. 

Fig. 5 The dependence of J{M) on Jc{M) for various values of T, where Jc{M) = J{M, T) 
at T = Tr. 



Fig. 6 Critical behavior of where a(T) is defined in —C, is the slope of the 

straight line. 

Fig. 7 The dependence of K{M) on M for various values of T. 

Fig. 8 The dependence of K{M) on K(.{M) for various values of T, where Kc{M) = 
K{M, T) at T = Tc. 



Fig. 9 The dependence of (3{T), defined in (pTl) , on T for various values of A. The slope 
of the straight line is —k = —2.2. 

Fig. 10 The dependence of the critical exponent n on A. 

Fig. 11 The direct relationship between /3(T) and a(T), when T is varied. 

Fig. 12 The direct relationship between K{M,T) and J{M,T), when M is varied but T 
fixed. 

Fig. 13 Same as Fig. 12 but in rescaled variables. 

Fig. 14 The probability distribution of the average density p in a configuration for three 
fixed T, generated by 10^ configurations at each T. The scale of p is arbitrary. 

Fig. 15 The dependence of J{M,p) on M for various values of p. 

Fig. 16 The dependence of J{M,p) on Jo(M) for various values of p, where Jo(M) = 
J{M,p) at p = pq. 

Fig. 17 Critical behavior of a{p), where a{p) is the slope of a straightline fit in Fig. 16. 
Fig. 18 The dependence of K{M, p) on M for various values of p. 

Fig. 19 The dependence of K{M,p) on Ko{M) for various values of p, where Ko{M) = 
K{M,p) at p = po. 
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Fig. 20 Critical behavior of P{p), where P{p) is the slope of a straightline fit in Fig. 19. 

Fig. 21 The direct relationship between f3{p) and a{p), when p is varied. 

Fig. 22 The effect of randomization on (a) Oi{T) and (b) where u is the number of 
steps of random walk. 

Fig. 23 The effect of randomization on (a) P{T) and (b) k, where v is the number of 
steps of random walk. 

Fig. 24 The effect of randomization on the relationship between P{T) and a{T). 
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